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ABSTRACT 



We detect dipping activity /modulations in the light curve of the four LMXBs in 
the 3-10 keV and 20-40 keV energy ranges. The spectral parameters derived from 
the fits to the INTEGRAL data are consistent with hot coronal structures in these 
systems where we find a range of plasma temperatures 3.0-224.9 keV. The unabsorbed 
X-ray to soft Gamma-ray fiux between 4-200 keV are 5.9x10"^'' erg s~^ cm~^ for 
XB 1916-053, 3.3x10-1° erg s'^ cm'^ for XB 1323-619, 21.6x10"^° erg s'^ cm'^ 
for X 1624-490 and 11.0x10"^° erg s'^ cm'^ for 4U 1746-371. The optical depth to 
Compton scattering, r, varies in a range 4.4-0.002 consistent with electron densities rig 
< 1.4x10^^ cm~^. In general, we find no significant difference in the dip and non-dip 
spectra in the ISGRI energy range (above 20 keV) for all the four sources. We only detect 
absorption differences between dipping and non-dipping intervals for XB 1916—053 and 
X 1624—490 in the JEM-X energy range. Fits in the 4-200 keV range including an 
additional photo-ionized absorber model for the two sources show that XB 1916—053 
has the highest ionized absorber amoung the two. 

Subject headings: X-ray: binaries - accretion, accretion disks - Radiation Mechanisms: 
general - X-rays: individual: XB 1916-053, XB 1323-619, X 1624-490, 4U 1746-371 



1. INTRODUCTION 

Some galactic low-mass X-ray binaries (LMXBs) show periodic dips in their X-ray intensity. 
The dipping behaviour repeats with the orbital period of the systems and is caused by periodic 
obscuration of the X-ray emitting region by a structure located in the outer regions of a disk particu- 



larly belived to be the impact region of the accretion flow from the companion star (jWhite &: Swank 



19821 ). The depth, duration and spectral properties of the dips vary according to the source and 
binary phase. There are about 10 members (some candidates) of these systems. The 1-10 keV 
spectra of most of the dip sources become harder during dipping. This change is accepted to be 



^Middle East Technical University, Physics Department, Inonu Bulvari, 356531, Ankara, Turkiye; 
solen@astroa.physics.metu.edu.tr 

■^Astrophysics Missions Division, Research and Scientific Support Department of ESA, ESTEC, Postbus 299, 
NL-2200 AG Noordwijk, The Netherlands 



- 2 - 



inconsistent with an increase in the photo-electric absorption by cool material since an excess of 
low-energy photons is also observed to be present in some of the systems (e.g., EXO 0748— 676: 
Cottam et allboOlh . 



Several approaches is used to model this spectral evoluti on (persistent-non-dip to dip spectra) 



The first approach is th e " absorbed plus unabsorbed" (e.g., iParmar et al.lll986l : ICourvoisier et al 



19861 : ISmale et al.lll992l ) where the persistent (non-dipping) spectral shape is used to model spectra 
from dipping intervals. The spectral evolution during dipping is accounted for by large increases 
in the column density of the absorbed component, and decreases of the normalization of the un- 
absorbed component. The latter decrease is attributed to electron scattering in the absorber. The 
changing normalization was difficult to justify and thus, a more complete and updated version of 
the model abov e was adopted as the "progressive cove ring", or by name "complex continuum" 
approach (e.g., IChurch et al.l 119971 : iBarnard et al.ll200ll ). where the X-ray emission is assumed to 



originate from a point-like blackbody (the neutron star), or disk-blackbody component, together 
with an extended power-law component. This approach uses a model where the extended compo- 
nent is partially and progressively covered in time by an opaque absorber to explain the spectral 
changes during dipping episodes. In such episodes, the blackbody component is also absorbed. 

More recently, the improved sensitivity and spectral resolution of XMM-Newton, Chandra and 
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20061 ) have recovered highly ionized Fe absorption features from these systems and demonstrated 
that modeling the changes in both the observed X-ray absorption lines and continuum seen during 
the dips and persistent emission can be explained by an increase in the column density and a 
decrease in the amount of ionization of a highly photoionized absorbing plasma. 

The dipping LMXB sources have been well stud i ed with the RXTE and Beppo-Sax missions in 



an energy range out to ^ 20 keV ( Parmar et al. 19991: iBaluciiiska-Church et al. 
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). Though both missions have harder energy range detectors (RXTE- 
HEXTE detector and Beppo-Sax-PDS detector), the diminishing sensitivity and short observation 
durations have not yielded effective analysis of light curves and spectra above 20 keV. This in the 
end yields spectral parameters that vary according to the detector and mission regardless of intrin- 
sic source variability. The ESA's 7-ray observatory INTEGRAL will allow for the determination 
of the widest energy band observed for dipping LMXB sources with the highest possible sensi- 
tivity above 20 keV. The energy range of the data (~ 15-500 keV) provided by the INTEGRAL 
observatory ISGRI detector supported by the lower energy JEM-X data should reveal the best 
parameters at high S/N that describe the broad-band spectrum of the dipping LMXBs together 
with the effects of scattering at high energies. Furthermore, we extend the investigation of the role 
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of photo-ionized plasmas via application of warmabsorber models in dipping LMXBs beyond the 
10 keV upper energy limit of pervious studies using data from INTEGRAL observatory. 

In this paper, we derive light curves and spectra for four dipping LMXB systems, XB 1916—053, 
XB 1323-619, 4U 1746-371, and X 1624-490. We discuss the dipping activity visible in the light 
curves and describe the broad-band time-averaged spectral characteristics of each source. We 
compare the dip and persistent (non-dip) spectra of the sources. We discuss and elaborate on the 
results of the spectral fits and compare these with previous findings in the light of the present 
physical models of accretion disk coronae and photo-ionized absorbers. 



1.1. Source properties 



The properties of the four dip sources studied here (XB 1916-053, XB 1323-619, 4U 1746-371 
and X 1624—490) are summarized in Table [TJ 



XB 1916—053 is the dip source with the shortest period known of 0.83 hr ([Walter et al.lll982l ). 
The presence of dips and the lack of X-ray eclipses fr om the companion star indicate tha t the 
system is viewed at an inclination angle, i, o f ~60-80° (IFrank et al.l 119871 : ISmale et al.lll988l ). By 



combining archival and R-XTE observations, IChou et al.l (|200ll ) were able to derive a precise dip 
recurrence interval of 0.83351411 (25) hr. The optical counterpart shows a modulation wi th a period 
1% longer (ICallanan et al.lll995l . and references therein). The presence of super humps (ISchmidtke 



1988 



Retter et al 



2002) have been proposed to account for this difference. XB 1916—053 has shown 



a 270 Hz highly coherent oscillation durin g a type I X-ray burs t indicating that the neutron star is 
spinning rapidly with a period of 3.7 ms dOallowav et allboOlh . The source has also shown quasi- 
perio dic oscillations (QPOs) at various frequencies ranging from ~0.2 to 1300 Hz (jBoirin et al 



200d ). The BeppoSAX 0.2-200 keV non-dip (or "persistent") spectrum may be modeled using an 
absorbed blackbody with a temperature, kT, of 1.62 it 0.05 keV and an absorbed cut-off power-law 
with a photori inde x, F, of 1.61 it 0.01 and a relatively high cut-off energy, Ecut, of 80 it 10 keV 
(jChurch et al.lll998l ). Narrow absorption f e ature s due to Fe XXV and FeXXVI in the XMM-Newton 
spectrum were discovered by lBoirin et al.l (120041 ) which was interpreted as the existence of a warm 
absorber in the system. On the other hand, for XB 1916—053 progressive covering gives a good 
description of the spectra l evolution in dipping from both ASCA and SAX observations, as well 
(|church et al.lll997l . Il998l ). 



XB 1323—619 exhibits X-ray bursts and irre gular intensity dips t hat repeat every 2.94 hr. The 

and the 



source was first detected by Uhuru and Ariel V (IForman et al 



1978; 



Warwick et al 



1981 



dips and bursts were discovered using EXOSAT (Ivan der Klis et al.l Il985l : iParmar et al.Ml989al ). 
During the dips, which typically last for 30% of the orbital cycle, the 1-8 keV intensity varies 
irregularly with a minimum of ~50% of the average value outside of the dips. The presence of 
periodic dips and the abs ence of X-ray eclipses indicates that the source is also viewed at i = 
60-80° (IFrank et al .1119871 ). The dip recurrence interval is not so precisely known, with the best 
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measurement being that of lBalucinska-Church et al.l (119991) of 2.94(2) hr. Th e source exhibits ~1 Hz 
QPOs which have been attributed to a disk instabihty (jjonker et al.lll999l ). The BeppoSAX 1.0- 
150 keV spectrum of the XB 1323—619 persistent emission can be modeled by a cutoff power-law 
with r = 1.48±0.01 and Er,,t = 44.1+^;j | keV t ogeth er with a blackbody with kT = 1.77±0.25 keV 



(jBalucihska-Church et al.lll999l ). iBoirin et al.l (j2005l ) examined the changes in the properties of the 
Fe XXV and Fe XXVI absorption features during persistent and dipping intervals. They found 
evidence for the presenc e of less-ionized mat erial in the line of sight during dips supporting the 
warm absorber scenario. I Church et al.l (|2005l ) also explains the XMM observations of this source 
using the progressive covering model where they attribute the absorption features to an existing 
accretion disk corona. 

4U 1746— 371 is located in the globular cluster NGC6441. The source exhibited 3 apparently 
ener gy independent inten sity dips separated by 5.0(5) hr, as well as flaring and bursting behav- 
ior (jParmar et al.l Il989bl ). The dips are shallow (25% of the continuum intensity) and largely 
energy in dependent in the usable energy range of the R-XTE Proportional Counter Array (2.1- 
16 keV) (IBalucihska-Church et al.ll2004l ). The best estimate of the dip period is 5.16 ± 0.01 hr 
(|Balucihska-Church et alll2004l ). The BeppoSAX 0.3-30 keV spectrum of 4U 1746-371 is unusu- 
ally soft and is dominated by a disk blackbody with kT = 2.82 =fc 0.04 keV, togeth er with a cutoff 
power-law with F = -0.32 ± 0.80 and Ecut of 0.90 ± 0.26 keV dParmar et al.lll999l ). 



X 1624—490 exhib its irregular dips in X-ray intensity that repeat every orbital cycle of 21 hr 
(| Jones fc Watsonlll989l ). This long orbital period means that as well as being the most luminous 
dip source, X 1624— 490 also has the largest stellar separation and is often referred to as the "Big 
Dipper". The X 1624— 490 X-ray light curve exhibits promine nt flares at energies >15 keV and 
dipping is evident at energies up to 6 keV dSmale et al.1 boOll ). These authors determine a dip 
recurrence interval of 20.8778(3) hr. In co ntrast to the other dip sources studied, no X-ray bursts 
have been observed from X 1624— 490 (e.g.. iBaluciiiska-Church et al.ll200d : ISmale et al.ll200ll ). The 
1-100 keV BeppoSAX persistent spectrum of X 1624—490 may be modelled using a blackbody 
with kT = 1.31 ± 0.07 keV an d a cutoff power-law with F = 2.0 and an ^^cut of ~12 keV 
(jBalucihska-Church et al.ll2000l ). Narrow absorption line feature s identified wi t h Fe XXV and Fe 
XXVI were discovered from X 1624-490 using XMM-Newton by IParmar et al.1 ([2002). 



2. Observations and data reduction 



The INTEGRAL (jWinkler et al.l 120031 ) pay load consists of two 7-ray instrum ents, one opti 



mized for 15 keV to 10 MeV high-resolution imaging (I BIS: lubertini et al.l I2OO3I ) and the other 



for 18 keV to 8 MeV high-resolution spectroscopy (SPI: IVedrenne et al.l 120031 ). IBIS provides an 
angular resolution of 12' full-width at half-maximum (FWHM), a fully coded field of view (FOV) 
of 8?3 X 8°, and a FWHM energy resolution, E/AE, of 12 at 100 keV. SPI provides an angular 
resolution of 2?5 FWHM, a fully coded hexagonal FOV of 13?2 x 13?2 and a FWHM E/AE of 
430 at 1.3 MeV. The extremely broad energy range of IBIS is covered by two separate detector 
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arrays, ISGRI (15- 1000 keV) and PI CsIT (0.175-10 MeV). In ad dition, the payload incl udes X-ray 
(JEM-X; 3-35 keV: lLund etallboosl l and optical (OMC; V-band: [Mas-Hesse et al]l2003 'l monitors. 
JEM-X has a fully coded FOV of 4?8 diameter and an angular resolution of 3' FWHM. All the 
instruments are co-aligned and operated simu ltaneously. The d ata used here were obtained from 
the INTEGRAL Science Data Centre (ISDC: [Courvoisier et al] [2003) archive. For this study, all 
publically available pointings between January 2003 and November 2004 (revolutions 30 to 250) 
where one of the selected sources was within 10° of the centre of the FOVs were included. Table [2] 
is an observation log. 

For this study, data from the ISGRI and JEM-X instruments were used due to their good 
sensitivity and broad energy coverage and processed using the Off-line Scientific Analysis (OSA 5.1) 
software provided by the ISDC. This includes pipelines for the reduction of INTEGRAL data from 
all four instruments. ISGRI and JEM-X use coded masks to provide imaging information. This 
means that photons from a source within the field of view (FOV) are distributed over the detector 
area in a pattern determined by the position of the source in the FOV. Source positions and 
intensities are determined by matching the observed distribution of counts with those produced 
by the mask modulation. After the standard data pipeline processing (dead time correction, good 
time-interval selection, gain correction, and energy reconstruction), source spectra and background 
subtracted light curves were derived using standard tools in OSA 5.1 . The positions of all the 
detected sources were fixed at the values given in the ISDC Referenc e Catalogue (version 23.0) 
which contains all known high-energy sources (see lEbisawa et al.l |2003| ) in order to derive reliable 
source fluxes. In addition, all sources that are at and above 6 sigma detection limit were extracted 
simultaneously during the analysis to avoid incorrect photon attribution between sources in order 
to map the field correctly. The coverage is less for JEM-X due to its smaller FOV compared to 
ISGRI. For further data manipulation we used standard software tools such as FTOOLS 6.0.4- 
6.1.1 and XSPEC 1 2. 2.1-12.3. . Pre liminary analysis of light curves and spectra can be found in 
Balman et aD ^20oi ): ISalmanl ^200^ ). 



3. Results 
3.1. Light Curves 

ISGRI and JEM-X light curves for all the four sources were derived using the standard OSA 5.1 
extraction algorithms. Event times were corrected to the solar system bary center when necessary. 
The data were grouped in bins of 10 s (XB 1916-053), 180 s (XB 1323-619), 120 s (X 1624-490) 
and 10 s (4U 1746—371) for the JEM-X analysis. The grouping time was increased to 240 s 
(XB 1916-053), 240 s (XB 1323-619), 240 s (X 1624-490) and 300 s (4U 1746-371) for the ISGRI 
analysis in order to increase the signal-to-noise ratio of the light curves. A search was made for 
any bursting activity since XB 1323—619, 4U 1746—371 and XB 1916—053 are known X-ray burst 
sources. For XB 1916—053, we detected 14 X-ray bursts with JEM-X, one of which is double 
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Fig. 1.— JEM-X (3-30 keV) and ISGRI (20-100 keV) folded light curves of XB 1916-053. The 
energy ranges are indicated on the panels. The folding period is given in Table [H 
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Fig. 2.— JEM-X (3-30 keV) and ISGRI (20-40 keV) folded light curves of XB 1323-619. The 
energy ranges are indicated on the panels. The folded light curve in the 40-100 keV energy range 
is not shown because of the low statistical quality. The folding period is given in Tabled) 



peaked. For XB 1323—619 some bursts with very low statistical quality were detected and removed 
from the data. No bursts were detected from 4U 1746-371 or X 1624-490 in either the JEM-X 
or ISGRI data. Some high intensity data points that appeared in X 1624—490 during non-dipping 
intervals have been removed while constructing spectra. These are most likely a result of few flaring 
activities of the source. For XB 1916—053 and X 1624—490 which have well determined dip periods, 
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Fig. 3.— JEM-X (3-10 keV) and ISGRI (20-40 keV) folded light curves of X 1624-490. The energy 
ranges are indicated on the panels. The folded light curve in the 40-100 keV energy range is not 
shown because of the low statistical quality. The folding period is given in Table [TJ 




Phose Phose 

Fig. 4.~ JEM-X (3-10 keV) and ISGRI (20-40 keV) folded hght curves of 4U 1746-371. The 
energy ranges are indicated on the panels. Some of the phases that are not covered well or show 
low statistical quality have been disregarded from the folded light curves. The folded light curves in 
the 10-30 keV and 40-100 keV energy ranges are not shown because of their low statistical quality. 
The folding period is given in Tabled) 

the uncertainties on the dip occurrence times over the observation duration are small compared 
with their orbital periods so the light curves from these sources were folded on the dip periods given 
in Table [TJ For XB 1323—619 and 4U 1746—371 the dip periods are not well determined, but the 
error that should be quadratically added from binning of the light curves to have enough signal to 
noise does not allow any period error determination better than presently known. We, therefore, 
also used the periods given in Table [1] for folding. In addition, the phases of the two light curves 
are arbitrary since we can not lock our epoch with the previous ephemerides due to the error in 
the periods for these two sources. 

Overall, we find dipping/modulation of light curves in all sources (not necessarily in all energy 
bands) folded on the periods in Tabled) The modulation depth normalized to the mean count rate 
can be defined as % modulation depth ~ lOOx (max rate - min rate)/(max rate + min rate)/2. 
XB 1916—053 indicate a modulation of 18% in the 3-10 keV energy range. Above 10 keV the JEM- 
X observations do not show evidence for any dipping activity for the source (see Fig. d])- With 
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ISGRI, the 20-40 keV reveals about 25% modulation depth. XB 1916—053 reveals only a small 
modulation of the light curve < 15% at the orbital period in the 40-100 keV energy range (see 
Fig. [H). For XB 1323-619, the JEM-X data show modulation of 40% in the 3-10 keV energy range 
for a small range of phases. It shows no modulation in the higher energies with JEM-X. The ISGRI, 
20-40 keV, data show no variation on the X-ray dip period, as well (see Fig. [2]). We note that this 
could be both because of a physical effect (e.g., precession of the disk) or due to large error in 
period determination. X 1624— 490 shows dipping/modulation in the JEM-X light curve of about 
25% (see Fig. [3]). The ISGRI light curve of the source has a modulation of about 40%. We detect 
some orbital modulation from 4U 1746—371 in the 3-10 keV (JEM-X) which is less than 10%. The 
ISGRI light curve shows dipping with about 45% modulation depth (see Fig. S]). A summary of 
modulation depths at certain energies for each source can be found in Table [3l 

We performed Fourier analysis of the time series obtained from the JEM-X and ISGRI data 
in order to recover the known frequency of orbital modulations, search for new periodic signals, 
and detect any red-noise features in t he low frequency band for XB 1916—053, X 1624—490, and 



4U 1746-371. The Scargle Algorithm (|Scarglelll982l : IPress et al.lll992l ) has been applied in the 3-10 



keV range to analyze the JEM-X time series data. Top panels of Figure[5]show Scargle periodograms 
(no power spectra are averaged ). The X-ray dip periods are noted on each panel of the figure (The 
periods are listed on Table [TJ The high peaks in these panels other than the noted periods are 
false features of high noise. The detection limit of a period at the 3a confidence le vel (99%) is a t 



a power level of 16.8 for X 1624-490 and 15.3 for XB 1916-053 and 4U 1746-371 (|Scarglelll982l '). 
Therefore, the signals (at the dip periods) are at 3a level for XB 1916—053 and X 1624—490, but 
the signal is insignificant for 4U 1746— 371and thus, the Scargle periodogram has been omitted. We 
note that the length of the time series used for both X 1624—490 and XB 1916—053 does not yield 
better errors than the published ones, thus we avoid comparison of any derived periodicity with 
the previous detections. However, Figure [5] confirms our folded light curves. In order to correct for 
the effects of windowing and sampling functions on power spectra, synthetic constant light curves 
were created and the actual power spectra were checked against the synthetic power spectra for 
false features. In addition, we did not perform prewhitenning or detrending of the light curves. 

We also performed power spectral analysis of the 20-40 keV ISGRI light curves in a simi- 
lar way to the JEM-X lightcurves using Scargle algorithm. However, the periodograms yielded 
no significant peaks in the power spectra of the three sources. In order to find the origin of the 
modulations detected in the ISGRI light curves we calcula ted stacked and averaged power spec- 



tra using standard Leahy normalization (jLeahv et al.lll983l ). We have averaged 108, 199, and 11 



power spectra for XB 1916-053, X 1624-490 and 4U 1746-371 with 2048, 1024, and 16384 fre- 
quency bins in each PSD, respectively. We detect a wide red-noise hump around/in the location 
of the orbital/dip frequencies in these sources over 3a (99% Confidence level) detection level and 
particularly a read noise peak at the location of the X-ray dip period of 4U 1746— 371 (see Fig. [5] 
lower /bottom panel). We do not detect a particular flattening of the red- noise in the 20-40 keV 
energy band in the frequency range of the orbital/dip periods which should be expected from disk 
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fed systems (jBelloni et al.ll2002l ). We think that the modulations detected above the 20 keV range 
in Figs. 1-4 are affected from these red noise humps, but the peak in 4U1746— 371is clear and it 
shows the highest percent modulation in the 20-40 keV range. 





4U 1746-371 (20-40 keV) 




Frequency (Hz) 



Fig. 5. — Top panels show the power spectra of XB 1916—053 and X 1624—490 between 3-10 keV 
derived using the Scargle algorithm (i.e., Scargle Periodograms) . Source names are given on each 
panel. The X-ray dip periods are labeled on the panels and 2H refers to the second harmonic of 
the period. See text for 3cr detection limits for XB 1916—053 and X 1624—490. Bottom panel is the 
average 20-40 keV power spectra of 4U 1746— 371. The expected (white) noise level is subtracted 
and the lowest level on the y axis corresponds to the 3fT detection limit for 4U 1746—371. 



3.2. Spectra 

The JEM-X and ISGRI data presented in Table [2] were used to derive a time-averaged broad- 
band spectrum of each of the four LMXBs studied in this paper. This data set provides a unique 
test that could only be performed with INTEGRAL since it is the only high-energy mission that 
frequently observes the galactic bulge region with large field of view instruments. During the time 
span of our observations, we do not detect any prominent high- or low-states of any of the LMXBs 
in our sample. 
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We first performed s pectral analys is on combined JEM-X and ISGRI spectra using XSPEC 
version 12.2.1 and 12.3.0 (|Arnaudlll996l ). A constant factor was included in the spectral fitting to 
allow for a normalization uncertainty between the tw o instruments . This factor was constrained 
to be within the expected range (i.e., 0.6-0.7) as (see iKirsch et al.l (|2005l )) as the parameter was 
thawed during the fitting pocedure. Spectral uncertainties are given at 90% confidence level (Ax^= 
2.71 for a single parameter) and upper limits at 95% confidence. We regrouped the JEM-X energy 
channels (1-256) by a factor of 3-4 to increase the signal-to-noise in the spectra. We used linearly 
sampled 40 channel response matrices for the ISGRI data where the data were similarly regrouped 
by a factor of 4-6. After the regrouping process, energy channels that have low signal-to-noise ratio 
have also been ignored. To account for systematic effects we added quadratically a 1% uncertainty 
to each spectral bin as recommended in the OSA 5.1 cookbook. We note that the addition of 
1-2% uncertainty for statistical concerns only reduces the Xv values, but does not affect the best fit 
parameters or parameter ranges. 



We fitted composite mod e ls to the broad-band spectra consisting of a disk blackbody (jMitsuda et al 



1984 



Makishima et al. 



1986 



Mitsuda et al 



or a blackbody, t ogether with either a cutoff 



power-law or CompTT model ( TitarchuD 1994 : Hua Sz Titarchuk 1995). Table H] shows the best-fit 



spectral parameters and their uncertainty ranges at 90% confidence level for each composite model. 
The disk blackbody or blackbody model was used to account for the emission from the disk and/or 
neutron star. The emission from an extended scattering/Comptonizing region was accounted for 
by the XSPEC model CompTT, a more physical model compared with the generic cut-off power-law 
model that only involves a simple power law decay with a high energy cut-off depending on the 
energy of the scattering electrons. During the fitting procedure, we also included the wabs XSPEC 
model to account for the neutr al hydrogen absorp t ion. These values were fixed at the values derived 
from the XMM-Newton data (iBoirin et al.ll2005l : iDiaz Trigo et al.ll2006l ) (also compared with the 
Nh values derived from the Beppo-SAX and R-XTE data to make sure all of them are similar). 
Figure [6] displays the fitted broad-band spectra of the four LMXBs. 

We also applied a single model to fit the four broad-band spectra. We used four models; the 
cut-off power law, CompTT, blackbody, or the disk blackbody for the analysis and found that none 
of the four LMXB spectra are consistent with a single emission model yielding reduced values 
>2.5 and/or nonphysical parameters except for XB 1323— 619. A fit of a power-law model with 



r3.0 



X 10- 



'^^ cm ^ with a power law 



absorption to the spectra of XB 1323—619 yields an Njj of 5.5_2 g 
photon index of 2.9lQQg and a normalisation of 0.6lg2^ photons keV~^cm~^s~^ at 1 keV {xl = 1-5). 
We also fitted the 4-200 keV spectra with a double blackbody or double disk-blackbody model. The 
fits resulted in reduced x^ values >2.0 and/or nonphysical parameters for all sources. 

The X-ray spectrum of X 1624—490 and XB 1323—619 as been modeled using a halo model 



because of the extensive dust in the viciiiity of the sources particularly of X 1624—490 (|Xiang et al 



20071 : ISmale et al.ll200ll : iDiaz Trigo et al.ll2006l : iBahicihska-Church et al.ll2000l : iBarnard et al.ll200ll ) . 

The radial intensity profiles of a source is expected to have excess over the point spread function 
due to the level of dust scattering at low energies mainly below 5-6 keV. We also applied the dust 
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model in XSPEC to the INTEGRAL spectra along with the two-component model of emission. 
This multiplicative model modifies a spectrum due to scattering off dust on the line-of-sight. The 
model assumes that the scattered flux goes into a uniform disk whose size has a 1/E dependence 
and whose total flux has a 1/E^ dependence. The spectral results in Table H] do not change with 
the inclusion of the multiplicative model. This is expected given the sensitivity, the large scale 
of the point spread function of INTEGRAL detectors and the energy range we are applying the 
model which is above 4 keV. We note that the detected blackbody temperatures (for the the two 
sources) are the same as the Beppo-Sax results including the dust contribution which means that 
the INTEGRAL spectra is not affected by dust scattering component. 

We also divided the JEM-X and ISGRI data into dipping and non-dipping episodes by selecting 
on phase intervals and created spectra for the LMXBs in our sample. The dip-phase bins in Figs. [T]- 
[3] (phases 0.9-1.05) were used together with the periods in Table [H The extracted dips are an 
average over different levels o f dippi n g, but the phase bins approximates the levels closest to the 



deepest dips iDiaz Trigo et al.l (j2006l ): iBarnard et al.l (j200ll ). We were able to recover statisticaly 



signiflcant and different dip and non-dip spectra for XB 1916—053 and X 1624—490. The analysis 
of the 4U 1746—371 data did not reveal any difference in the dip and non-dip spectra. Given the 
low count rates and the short dip intervals observed from XB 1323—619, we could not derive a 
statistically signiflcant dip spectra with JEM-X in 15 ksec (dipping time in the given exposure 
time) and the high energies (ISGRI) show no variation, hence no dipping activity at the given 
period (see Fig. [2|). 

Fig. [7] shows the simultaneously fit dip and non-dip spectra of XB 1916—053 and X 1624—490 
with only a composite model of a blackbody and compTT. The top right hand panel of Fig. [7] shows 
the fitted XB 1916—053 spectra and the top left panel shows the X 1624—490 spectra. The composite 
spectral model was simultaneously fit to the dip and non-dip spectra of the two sources for a better 
determination of the differences in the fitted results (the fitted model was wabs*(bbody+compTT)). 
The residuals of the top panels of Fig. [7] indicate distinctive deviations from the mean and the 
reduced values are >2.4 in the top panels (the d.o.f. are similar to the ones reported on 
Table H]). Acceptable fits were achieved when two different neutral hydrogen column densities for 
dip and non-dip spectra were assumed along with an additional absorption line at 6.9 keV for 
XB 1916—053 (dip spectra) and an emission line for X 1624—490 (non-dip spectra) in the fitting 
proccess (additional models of gabs and gauss were used for modeling the lines). The spectral 
parameters for the composite blackbody and compTT model of the dip and non-dip spectra for 
XB 1916—053 and X 1624—490 resemble the best fit values on Table H] very closely, so we do not 



inclu de them in this paragraph. On the other hand, we found a previously detected (jBoirin et al 
20041 ) absorption line in XB 1916-053 at 6.9151 ^eV. Addition of the hue improves the fluctuation 



of the (flt) residulas from above 2a down to less values in the vicinity of 6-7 keV. We also found 

hO.2 
-0.2 



a previously detected emission line in the non-dip spectrum of X 1624—490 at 6.8]']g 2 keV with an 



integrated flux of 2.1^*^0 gX 10 ^ photons s ^ cm ^. Addition of the line improves the fluctuation 
of the (flt) residulas from 3a down to less values in the vicinity of 6-7 keV. 
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The detection of an increased amount of absorption in dip spectra for both sources together 
with the presence of an absorption hne in XB 1916—053 suggests the existence of photoionized 
absorbers solely from the INTEGRAL data. However, INTEGRAL data, alone, does not account 
for existence of photoionized absorbers in the persistent emission. Fig. [7] middle and bottom panels 
show the non-dip and dip spectral fits to both XB 1916—053 and X 1624—490. In order to compare 
with the results from the XMM-Newton analysis of these two sources with models of photoionized 
absorbers, we applied the warm absorber model, warmabs (an additional model in XSPEC), com- 
piled with XSPEC 12.3.0 to both spectra of XB 1916-053 and X 1624-490. Warmabs models the 
absorption from a photoionized plasma in the line of sight. Column densities of ions (including 
small cross-sections) are coupled through a photoi onization model. Warmabs calculates spectra us- 



ing s t ored l evel populations calculated by XSTAR (IKallman fc Bautistall200ll : lBautista fc Kallman 



2001 



2000) assuming a given continuum spectrum. For the ionizing flux a photon index of F = 2 



was assumed in this study consistent with the X-ray observations of these systems. The spectral 
parameters of warmabs are an absorption column density of hydrogen of the photoionized plasma, 
the logarithm of the ionization parameter log(^), and the turbulent velocity broadening along 
with elemental abundances. The ionization paramter ^ = L/nR'^ where L is the luminosity of the 
ionizing source, n is the gas density, and R is the distance from the ionizing source. The middle 
panels in Fig. [7] show the non-dip and the bottom panels show the dip spectra fitted with the 
composite model of a neutral hydrogen absorption, a blackbody and a compTT including a warmabs 
model to account for the excessive absorption and/or the absorption line. The fits performed with 
the dip and non-dip spectra of XB 1916—053 and X 1624—490 yield reduced values around 1.0- 
0.83 (for 36-31 d.o.f.) . Using the fits to the INTEGRAL data, we derive a log(0 of S.Stoi for the 
persistent and 2.8l;g 2 ^^P spectra of XB 1916—053. The absorption column of the ionized 

absorber is 4.8]']}'3X 10^^ cm~^ during the dipping episodes and l.T^J'gX lO^'^ cm~^ during non- 
dipping intervals. The fits performed with the dip and non-dip spectra of X 1624—490 yield a log(^) 

I (\ cy _i_n o 

of 3.4_q'2 for the persistent and 2.6_g'2 for the dip spectra. The fit results indicate an absorption 
column of 4.0lj^ gX 10^^ cm~^ for the ionized absorber in the persistent and 15.8l^ gX 10^^ cm~^ in 
the dipping spectra. The parameter to the line of sight are kept constant as in Table H] during 
analysis. The pa.ramte r was fixed at the values determined from XMM-Newton observations 
(jPiaz Trigo et al.l 120061 ): fixed at the values of 100 km s^^ (dip) and 300 km s^^ (non-dip) for 
X 1624-490 and 200 km s-^ (dip) and 2000 km s^^ (non-dip) for XB 1916-053. 



4. Discussion 

We report the results of archival INTEGRAL observations of four dipping LMXBs, XB 1916—053, 
XB 1323-619, 4U 1746-371 and X 1624-490. The data presented in Table [2] were used to derive 
time-averaged broad-band spectra of the four LMXBs (January 2003-November 2004) which is a 
unique capability of the INTEGRAL observatory as the data analysis is performed by averaging 
observation science windows over long periods of time. We, also, investigated the dip and non-dip 
spectra. During the time span of our data, we do not detect a prominent high or low state of any 
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of the LMXBs in our sample. 



The LMXB spectra have been assume d using several different approaches. One of the ap- 
proaches, Western model (jWhite et al.lll988l ). is a two-component emission model where the harder 
X-ray component is the Comptonized emission originating from an ADC in the LMXB spectra and 
it has the form of a power-law with an energy cut-off at the limit of the Comptonizing electrons and 
only the very bright sources required an additional blackbody component from the boundary layer 
between the accretion disk and the surface of the neutron star. Another model, Eastern model 



(jMitsuda et al.lll989l ). assumes that photons from the neutron star are Comptonized in a small cen- 
tral region on the disk and harder component orig inates from a multi-temperature dis k blackbody 
model. A descendent model, Birmingham model ([Church &: Balucinska-Churchlll995l ). takes after 
the Western model, but assumes neutron star blackbody for the soft components. High-resolution 
spectroscopic observations of LMXBs have shown in some cases the presence of extended X-ray 
emitting plasmas or onflows surrounding the accretion disk. The detection of these components 
and their kinematic properties has been made possible by the observation of emission or absorp- 
tion lines from very highly ionized metals. This also diagnoses several different e mitting zones in 
LMXBs such as disk atmospheres and coronae ()Jimenez-Garate et al.l 12002 . |2003| ). In this paper, 
we fitted spectra assuming several different models and constructed Tabled] to test the INTEGRAL 
data given the presumed approaches to derive the spectral parameters, also, introducing a warm 
absorber model to model the absorption over 4-200 keV in the light of recent model construction 
for dipping LMXBs. 

This INTEGRAL study uses better S/N spectra and lightcurves incomparison with other stud- 
ies, particularly above 20 keV for the four LMXBs in question. None of the RXTE studies of these 
four sources include usable spectra or light curve above 20 keV. The Beppo-Sax has the opportu- 
nity to extend out to 100 keV using the PDS instrument which has been used in the literature. 
However, the shorter exposures 15-30 ksec and low count rates, do not allow for a high S/N spectra 
above 20 keV, the effective area also falls readily and there are no PDS light curves available in the 
literature above 20 keV. In general, for two of the sources 4U 1746—371 and X 1624—490, there is no 
effective study of the spectra or light curve above 15-20 keV. For XB 1916—053 and XB 1323—619 
the PDS spectra have been helpful, however the S/N is not favorable to the INTEGRAL study 
here. INTEGRAL countrates are three to four times better for above 500 ksec integration time (we 
used about 1 Msec exposure on each source with ISGRI). As a result ISGRI results in this study 
provide new information on these objects and they are supported by the JEM-X detector (to avoid 
cross-calibration issues that could arise from using other observatory data) to provide large energy 
band coverage. 

Our spectral results in Table |4] yield blackbody temperatures of 1.7 keV and 1.4 keV, for 
XB 1323—619 and X 1624— 490 which are in excellent agreement with the BeppoSax results (seed]) 
within statistical errors. On the other hand, XB 1916—053 and 4U 1746—371 yields different black- 
body temperatures of 2.3 keV and 0.6 keV, respectively where the BeppoSax results reveal 1.62 
keV and 0.5 keV for the same sources (see [1]). The INTEGRAL result and BeppoSax results 
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for 4U 1746—371 agree within statistical errors. Moreover, the Wien photon temperature kTo 
(parameter-seed in photons) varies in a range 0.4-2.5 keV in the four sources (even if we derive 
similar blackbody temperatures). These differences could be attributed to the existence of more 
than one blackbody component for at least the cases like XB 1916— 053 and 4U 1746— 371. At 
suffciently high mass accretion rates, the boundary layer and the accretion disk are both radi- 
ating in the optically thick regime with kT 1-2 keV. It is expected on theoretical grounds and 
demonstrated that the boundary layer has higher temp erature than the accretion disk yielding two 
different blackbody temperatures in the X-ray spectra (jGilfanov et al.l 120031 : iGogiis et al.ll2007l ). 



We find that the normalisation parameter for the disk blackbody model, proportional to 
Rf^cosO, yields fairly small real inner disk radii of less than 14.5 km, 5.2 km, and 12.8 km (all 
at 3a maximum hmit) for XB 1323-619, XB 1916-053, and X 1624-490, respectively {9 is the 
inclination angle suitably taken between 60°-80°). We stress that at 2a limit these radii are less 
than a neutron star size. On the other hand, this value is between 27-37 km for 4U 1746—371. As a 
result, a disk blackbody emission is very unlikely for the three sources except for 4U 1746—371. We 
note t hat the 3a maximum limits on apparent inner disk radii are corrected to real radii accord- 
ing to iKubota et al.l (|l998l ) . We caution that normalisation of the disk blackbody model involve 
distance in the calculation which results in a slight uncertainity. Radii change only linearly with 
distance and only a factor of 2-3 times (in kpc) difference in distance will effect the general trend 
discussed above. In comparison with this, we derived the size in radii of the blackbody emitting 
regions for the four sources. We calculate that for XB 1323—619 and XB 1916—053 the radii is 
between 0.8-1.02 km. For X 1624— 490, this radius is about 9.7-11.1 km (e.g., the size of a neutron 
star). 4U 1746— 371 shows an emitting region of 37-51 km which is physically inconsistent with 
emission from the surface of a neutron star. These ranges above correspond to 3a error ranges of 
the used parameters. INTEGRAL results suggets that XB 1323—619 and XB 1916—053 accrete to 
a small region on the surface of the neutron star and 4U 1746—371 shows only a disk blackbody or 
boundary layer emission. The disk blackbody temperature of 4U 1746— 371detected in this study 
is consistent with the BeppoSax results as well (see [1]) with only slight difference in disk blackbody 
temperature. It is difficult to make such deductions for X 1624—490 using solely INTEGRAL data 
except that a disk blackbody emission is inconsistent with the source. In addition, we note that the 
normalisation of the blackbody model involve distance in the calculation resulting in some uncer- 
tainity. Radii change only linearly with distance. As a result, a factor of only 10 times difference 
in distance for the case of XB 1916—053 and XB 1323—619 and a factor of only 2-3 difference in 
the case of X 1624—490 will effect the general trend discussed above. 

A comparison with the spectral results obtained from the BeppoSax and R-XTE observations 
using the Cut-off power-law model for the second emission component shows that the time average 
data yield somewhat different results for t he E„,t and the pho ton index parameters (see Table [H 
for individual r esults). For XB 191 6-053, Ichurch et aP (|l998l ^ find F = 1.61±0.01 and E^^t = 



80ibl0 keV and lNarita et al.l (j2003bl ) find a range of powerlaw indicies between 1.74-2.05 and Ecut 
of > 30 keV (90% confidence limits). Our results are consistent with R-XTE observation analyses. 
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For X 1624-490, Balucinska-Church et alJ ((20001) find T = 2.0±0.06 and E^ut >12 keV where we 
find a photon index > 0.84 a r id a lo w ii^cut value of 1.23-4.6 keV including uncertainty ranges. For 
4U 1746-371, IParmar et ail ijlddd ) find F < 0.5 and Ecut = 0.90±0.26 keV where we find F = 



2.18±0.8 and E'cut value of 6.6-11.2 keV. For XB 1323-619, iBalucihska-Church et al.l (|l999l l find 
F = 1.48±0.01 and ^;cut = 44±5.1 keV (BeppoSax results) and Ba rnard et aU (|200lh derive F = 
1.23±0.07 and E^nt >44 keV f R-XTE results) . Our results yield F = 3.0±0.2 and an unbound cut- 
results show an F = 1.89-1.91 for Ecut > 117, 



off energy larger than 156 keV. lBarnard et al. 



(2001 



however without inclusion of a blackbody model (90% confidence limits). These are discrepancies 
regardless of error estimates and such differences between the -Ecut and the photon index parameters 
exists also between BeppoSax and R-XTE results. Apart from the fact that there are possible cross- 
calibration inconsistencies, and that the sources could be variable, this could also be attributed to 
the Cut-off power-law model being an inadequate model of emission to represent the data in the 
harder energies particularly in the ISGRI energy band (e.g., relativistic effects). This strengthens 
the fact that the Comptonizing plasmas in these systems needs to be modeled more properly. 

We applied the CompTT model in XSPEC inorder to account for the Compton scattering regions 
in these systems and found a range of plasma temperatures of 18.8lg'{j^ keV for XB 1916—053, 
196.4128-5 ]jgv for XB 1323-619, 9.2^1^ keV for X 1624-490 and 4.9t|^ keV for 4U 1746-371 given 
the best fit results of time-averaged data. These temperatures indicate persistent hot Comptonizing 
regions, accretion disk coronae (ADCs), in these systems. The optical depth of Compton scattering 



regions can, also, be derived from the fits, r, is in a range 0.9ll^Qg 



for XB 1916-053, 0.006 



XB 1323-619, 0.2+^;2 for X 1624-490 and l.i:l:o;9 for 4U 1746-371 which, in addition, indicates 
that the Comptonizing regions are relatively tenuous. XB 1323—619 has the low est value. Using 
the model parameter d escriptions of CompTT, Ue^r^/radfT at high optical depth (lTitarchuklll994l : 
Hua &: Titarchuklll995l ) , we derived maximum limits for the electron densities in the X-ray emitting 
regions ofthe four systems as: XB 1916-053, ne<7.6xl0^'^ cm"^, 4U 1746-371, ne<1.4xl0^^ cm'^ 
X1624-490, ne<l. 1x10^5 cm'^, and XB 1323-619, np<5.4xl0^ cm"^. W e have assumed 0.25 xr^d 
for the emitting region (r^d taken from lChurch &: Balucihska-Churchll2004l ) . If the emitting region is 
O.lxr^d the densities increase by a factor of four. If it is spread over the entire disk, the densities are 
a factor of four less. The electron densities that can be calculated fr om the ionizatioii param eters 
derived from the fits with photoionized absorbers for these systems in iDiaz Trigo et al.l (|2006l ) yield 



2.3 



-0.001 
-0.004 



for 



a range of 8-0.2 x 10 cm^ for the ionized absorbers in these systems. These values are consistent 
with the maximum limits in the Comptonizing regions. Then, the suggested warm absorbing regions 
may have similar densities with the existing coronal regions. 

The Compton radius in an accretion disk corona (ADC) with temperatur e kT is rc— GMmp/kTp^asma 
(cf. , IChurch &: Balucihska-Churchll2004l : iBalucihska-Church &: Churchll2005l) where the hydrostatic 



equilib rium fails. Outside this radius, the corona will dissipate as a wind. IChurch &: Balucihska-Church 



(|2004l ) shows that the Compton radii for some LMXB dippers and the ADC radii calculated from 
the ingress and egress timescales in their lightcurves are very similar. We calculated the Compton 
radii for the four LMXBs given the plasma temperatures from the fits with the CompTT model as 
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follows : XB 1916-053 has rc=il.3xl0^° cm, 4U 1746-371 has rc~1.8xl0^° cm, X 1624-490 has 
rc~3.9xl0^*^ cm, and XB 1323—619 has rc~l-lxlO^ cm. The extent of the Comptonizing regions, 
except for XB 1323—619, are large and it could cover 13-77% of the disk (by area) given the statis- 
tical errors on the plasma temperature. For XB 1323—619 the coverage of the ADC is only about 
1% of the disk (by area), most likely a small region confined to inner parts of the disk (given the 
size of calculated rc). 

We have also modeled the absorption in the broad-band INTEGRAL energy range with pho- 
toionized absorbers owing to the detection of mainly iron absorption lines from these systems 
(other absorption lines have also been detected). We were able to satisfactorily separate dip and 
non-dip spectra from two of the systems and studied, for the first time, 4-200 keV spectra modeling 
photoionized absorbers in these systems. A comparison of the spectral results with the XMM- 
Newton analysis in the 0.2-15.0 keV band (see Table yields log(^) values larger in the case of 
XB 1916—053 and slightly smaller in the case of X 1624—490 for both non-dip spectra obtained from 
the time-averaged INTEGRAL data. For dip spectra both sources show slightly higher ionization 
parameter (compared with deepest dips in the XMM-Newton analysis). This shows that INTE- 
GRAL detects a more highly ionized absorber in XB 1916—053 (even higher than X 1624—490) 
which was not recovered in the 0.2-15.0 keV analysis. The neutral absorption of the highly ionized 
absorber is significantly lower compared with the XMM-Newton results at all times. The XMM- 
Newton analysis, also, assumes a photon index of F = 2 for the ionizing spectrum as assumed in 
this analysis, but with an additional cut-off energy of 86 keV and 12 keV for XB 1916—053 and 
X 1624— 490, respectively. A comparison with the Chandra (0.3-10.0 keV) spectral results of the 
non-dip emis sion shows that INTEGRAL detects the same log(^) parameter. A value of 3.73±0.13 



laria et al 



is derived bv lJuett &: Chakrabartyi ([200a) assuming a power law of -1 for the incident ionizing flux. 



For X 1624-490 
which is higher then the INTREGRAL result. 



(|2007l ) measures approximately a log(^) of 4.1 for the non-dip spectra 



5. Summary and Conclusions 

We have investigated 4-200 keV broad-band spectrum and temporal characteristics of four 
dipping LMXB systems, XB 1916-053, XB 1323-619, X 1624-490, and 4U 1746-371. INTEGRAL 
has shown us for the flrst time the time-averaged spectral parameters of these systems together 
with the dip and non-dip spectra. This study also demonstrates for the first time, a high sensitivity 
study of spectra and timing characteristics of these four LMXBs above 20 keV (simultaneously with 
3-20 keV band data). INTEGRAL reveals the spectral characteristics of the hottest regions in these 
dipping LMXBs together with the nature of the absorbers within the systems. We find that the 
persistent emission spectral parameters of X 1624—490 and XB 1916—053 remain mostly the same 
in and out the dipping intervals given the 4-200 keV energy range except for the changing absorption 
below 10 keV for the two systems. 4U 1746—371 shows a spectrum that is consistently the same, in 
and out the X-ray dips between 4-200 keV given the INTEGRAL sensitivity. INTEGRAL spectra 
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are consistent with existence of persistent hot Comptonizing coronal regions (ADCs), with plasma 
temperatures using Comptonized plasma models as 18.81^^ keV for XB 1916-053, 196.41^8-5 keV 
for XB 1323-619, 9.2^1^ keV for X 1624-490 and 4.915:9 ^eV for 4U 1746-371. We caution that 
the ADC temperature for XB 1323—619 is very high incomparison to the rest and even within the 
LMXBs and the tem perature limit of the CompTT model to satisfactorily represent the ADC falls 



readily after 50 keV (jFarinelli et al.l 120081 ). Using our fits, we derive a range of optical depth to 
Compton scattering r in these regions as O-dtol for XB 1916-053, 0.006l[j;|][}^ for XB 1323-619, 
0.22l^JJ for X 1624-490 and l.lloJ for 4U 1746-371. Given these results XB 1323-619 has, 
by a large factor, the hottest ADC with the lowest optical depth. None of the sources exhibit 
spectra consistent with a single model of emission except for XB 1323— 619(i.e., a single power 
law model can be used to model the spectra). Therefore, we stress that the spectra of these 
sources should never be fitted using a single model of emisison regardless of the energy band (i.e., 
as in Chandra or XMM-Newton bands). In addition, the four systems can not be fitted, alone, 
with a double blackbody emission or double disk-blackbody emission models, either. Fits to the 
INTEGRAL spectra suggests a third component particularly of blackbody nature (due to the the 
detected different seed-in photon temparature compared with blackbody temperature) that could 
be existent in the spectra. 

For the first time, we detect ~15%, 40%, and 45% modulation depth (variation normalized to 
mean flux) in the 20-40 keV energy range at the X-ray dip periods from XB 1916—053, X 1624—490, 
and 4U 1746— 371, respectively. We belive this could be as a result of scattering off of the bulge 
at the accretion impact zone (fixed in the orbital plane). Photoelectric absorption by material at 
the impact zone is not as desirable an explanation since after 20 keV Comptonization and energy 
independent scattering is expected to dominate. This is also supported by the large red noise 
(humps) found in the average power spectra above 20 keV. 

We have studied the effects of photionized absorbers in dipping LMXBs for the first time 
in a broad energy band 4-200 keV from XB 1916-053 and X 1624-490. We find that the most 
highly ionized absorber is in XB 1916—053 with an ionization parameter log(^) of 3.81q4 although 
X 1624—490 has more luminosity. This indicates the importance of system geometry and size. 
Moreover, inclusion of spectra from a harder energy band compared with the XMM-Newton or 
Chandra data help to constrain the absorber properties better. 

The INTEGRAL spectrum reveals a small size accretion region that is less than 1 % of the 
surface area of a typical neutron star in an LMXB derived from the blackbody fits for XB 1916—053 
and XB 1323—619. This is difficult to explain, alone, by the large amount of scattering in these 
systems. On the contrary, accretion region sizes either about the size of a neutron star or 16-20 times 
larger ( than a neutron star in terms of accretion surface area) are consistent with X 1624—490 and 
4U 1746— 371, respectively. Such results from INTEGRAL reveals different accretion geometries 
and emission regions for these dipping LMXB systems in contrast with the very similar modeling 
of such regions used in previous studies (i.e., XMM-Newton, RXTE, Beppo-Sax). 
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The INTEGRAL broad-band hard X-ray and 7-ray emission from the four dipping LMXB sys- 
tems show that these systems have different Comptonizing regions in size and in plasma temperature 
and has different accretion geometries and emitting regions. We calculate that Comptonizing coro- 
nal regions cover a larger fraction of the disk and are colder in comparison with the conventional 
ADCs in three systems (i.e., XB 1916-053, X 1624-490, 4U 1746-371), such regions vary with 
respect to the source and are not only strongly dependent on the central source luminosity. INTE- 
GRAL detects that ADC in XB 1323—619 is much smaller in size (much more hotter) in comarison 
to the other three systems and probably confined to the inner regions of the disk. We do not find a 
direct correlation of source luminosity with plasma temperature in the ADCs of the four systems. 
These extended highly ion ized coronal regions (ADCs) may act as an absorber to the incident radi- 
ation ([Church et al.ll2005l ). As a result, we would be viewing these systems through highly ionized 
absorbing regions at all times. On the other hand, INTEGRAL results confirms the existence of 
warm absorbing regions and that the absorbing material looks less photoionized during the main 
dips and more ionized as viewed in the persistent emission intervals in a broad energy band of 4-200 
keV for at least two of the dippers XB 1916—053 and X 1624—490. These absorbing regions can 
be associated with an accretion disk atmosphere (e.g., EXO 0748—676: Ijimenez-Garate et al.ll2003l ) 
However, the two systems are not well known X-ray line emitters. It is unclear, at this stage, if the 
absorber is the same in a given location, or if different absorbers are viewed during dipping and 
persistent emission episodes. It may be that it is a combined effect of different absorbing regions 
thoughout the disk (e.g., inner disk and outer disk) superimposed on the spectra. 
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Table 1: Properties of the studied LMXBs. L^r is t he 0.6-10 keV luminosity in units of 10^^ erg s 



for distances d ta ken fromlDiaz Trigo et alJ (120061 ). except for XB 1323—619 where the 0.5-10 keV 



luminosity is from lBoirin et al. ( 2005 ). N^^'^ is the persistent and dipping emission column density 
in units of 10^^ atom cm~^ of the best-fit photo-ionized abs orber with an ion i zation parameter 



log(£ ) (in units of erg cm s ^) measured by XMM-Newton (jBoirin et al. 



2005; 



Diaz Trigo et al 



2p06i). Note that no such absorber was required with the fits to 4U 1746— 371 ( 
2006 ). 



Diaz Trigo et al 



LMXB 


(hr) 


^'36 
(erg s^i) 


d 

(kpc) 


Non-dip / Dip 


log(«) 
Non-dip / Dip 


Pdip 
Reference 


XB 1916-053 
XB 1323-619 
4U 1746-371 
X 1624-490 


0.833351411(25) 
2.94(2) 
5.16(1) 
20.8778(3) 


4.4 
5.2 
10.1 
47.5 


9.3 
10 
10.7 
15 


4.2l»j/28l^ 
3.81^4 / 14t^0 

13tL / 29tl 


o nt; + 0'04 ,r, c:c: + 0.04 
■^■"•^-OMI ^■•^■^+0.04, 
n Q+0.04 / o n+0-2 
■^■^-0.04 / •^•"-0.2 


Chou et al. (2001) 
Baiuciriska-Church et al. (1999) 
Baiucinska-Church et al. (2004) 
Smale et al. (2001) 



Table 2: Observation Log. The longer exposure times for ISGRI compared with JEM-X are due to 
the larger FOV of this instrument. Mean count rates are given in the energy range 4.0-20 keV for 
JEM-X and 20-200 keV for ISGRI. 



LMXB 




Exposure (ks) 


Count Rate (s ^) 


Tstart-Tstop (MJD) 






JEM-X 


ISGRI 


JEM-X 


ISGRI 


JEM-X 


ISGRI 


XB1916 


-053 


293 


915 


3.5 


1.2 


52952.4-52963.5 


52708.7-53310.6 


XB 1323 


-619 


497 


834 


0.96 


0.66 


53085.7-53263.9 


52650.7-53263.9 


4U 1746- 


-371 


316 


626 


6.4 


0.43 


52713.4-53068.1 


52713.3-53292.7 


X 1624- 


490 


398 


756 


9.4 


0.42 


52701.8-53062.0 


52700.7-53265.2 



Table 3: A summary of modulation depths for the dipping LMXBs analyzed in this work. The 
depths are given in percentage modulation normalized to mean count rate (~ 100 x (max rate - 
min rate) /(max rate -|- min rate)/2). 



LMXB 3-10 keV 20-40 keV 40-100 keV 



XB 1916-053 
XB 1323-619 
4U 1746-371 
X 1624-490 



18% 
40% 
<10% 
25% 



25% 
<14% 
45% 



<15% 




1 00 200 
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0.01 

10"- 
10"*^ 

10"; 



X 1524-490 -, 





Energy (keV) 





4U 1746-371 



Energy (keV) 



20 50 
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Fig. 6.— Simultaneous 4-200 keV fits to the JEM-X and ISGRI spectra of the four dipping LMXBs 
(4-20 keV; JEM-X and 20-200 keV; ISGRI). The crosses indicate the data, the solid lines show 
the summed model and the dotted lines show the individual models. The second panel under each 
count rate spectra display the residuals in standard deviations. Next to each count rate spectrum an 
unfolded spectrum in E^F is plotted to show the contribution of each component of the spectrum. 
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Table 4: Best-fit spectral parameters obtained using a cutoff power-law or compTT model together 
with either a disk blackbody or blackbody component. A^h is the absorbing column in units of 
10^^ atom cm~^, kT is the blackbody or disk blackbody temperature in keV, K-q and are the 
normalisations of the disk blackbody and blackbody models, F is the power law photon index, Ecnt 
is the cutoff energy in keV, Kcut is the model normalisation photon cm^^ s^^ keV^^ at 1 keV, F is 
the 4-200 keV flux in units of 10"^'^ erg cm~^ and L is the 4-200 keV source luminosities at the 
distances given in Tabled] in units of 10^^ erg kT^ is the input soft photon (Wien) temperature 
(keV), A;Tp is the plasma temperature (keV) and r is the plasma optical depth (electron scattering 
opacity). 
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Fig. 7. — Simultaneous 4-200 keV fits to the JEM-X and ISGRI time-average spectra of 
XB 1916-053 (left) and X 1624-490 (right) (4-20 keV; JEM-X and 20-200 keV; ISGRI). Crosses 
indicate the data, the solid lines show the summed model and the dotted lines show the individ- 
ual models. The second panel under each count rate spectra display the residuals in standard 
deviations. The top panels are dip and non-dip spectra fitted simultaneously with only the black- 
body and compTT model (wabs*(bbody+compTT)). The middle and bottom panels are non-dip 
and dip spectra, repectively, fitted with the composite model including the warm absorber model 
(wabs*warmabs* (bbody+compTT) ). 



